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Abstract — A planar quantum-well device made of a gapless graphene nanoribbon with edges in contact 
with gapped graphene sheets is examined. The size-quantization spectrum of charge carriers in an asym- 
metric quantum well is shown to exhibit a pseudospin splitting. Interface states of a new type arise from 
the crossing of dispersion curves of gapless and gapped graphene materials. The exciton spectrum is calcu- 
lated for a planar graphene quantum well. The effect of an external electric field on the exciton spectrum 
is analyzed. 
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1. INTRODUCTION 

The creation of graphene, a monolayer of carbon 
atoms forming a regular hexagonal lattice [1-3], has 
stimulated extensive experimental and theoretical stud- 
ies along various lines of research. Graphene' s unique 
properties make it a promising material for a new gen- 
eration of carbon-based nanoelectronic devices. In par- 
ticular, carrier mobility in graphene amounts to 2 x 10 5 
cm 2 /V-s, and ballistic transport is possible on a sub- 
micrometer scale [4, 5]. 

Over the past seven years, numerous theoretical 
and experimental results have been reported on elec- 
tronic properties of nanometer-wide ribbons of graphene 
(nanoribbons). Among the first were studies of elec- 
tronic states of graphene nanoribbons using the Dirac 
equation under appropriate boundary conditions [6, 7]. 
The electronic properties of a graphene nanoribbon 
strongly depend on its size and edge geometry [8]. In 
terms of transport properties, graphene nanoribbons 
are highly reminiscent of carbon nanotubes [9, 10] since 
free carrier motion inside them is also one-dimensional. 

A field-effect transistor (FET) based on a 2 nm 
wide and 236 nm long graphene nanoribbon was fabri- 
cated in a recent study [11] (nanoribbons of widths be- 
tween 10 and 60 nm were also studied). The graphene 
nanoribbon was made narrow enough to open a gap of 
width required for room-temperature transistor oper- 
ation. However, it is less compact than the graphene 
quantum-dot transistor 30 nm in diameter discussed in 
[12]. 

In this study, we examine a planar quantum-well 
device made of a graphene nanoribbon whose edges 
are in contact with gapped graphene sheets. 



A bandgap opening in graphene can be induced by 
several methods. First, graphene can be deposited on 
a hexagonal boron nitride (h-BN) substrate instead 
of a silicon-oxide one. This makes its two triangu- 
lar sublattices nonequivalent, inducing in a bandgap 
of 53 meV [13]. Second, epitaxial graphene grown on 
a silicon-carbide substrate also has a nonzero bandgap 
[14]. According to angle-resolved photoemission data, 
a bandgap of 0.26 eV is produced by this method [15]. 
Third, a hydrogenated derivative of graphene synthe- 
sized recently, graphane [16], has been predicted to 
have a direct bandgap of 5.4 eV at the T point [17]. 
Fourth, ab initio calculations have shown that Cr03 
adsorption on graphene induces a gap of 0.12 eV [18]. 
In the first two methods, a heterogeneous substrate can 
be used, such as an h-BN — Si02 nanoribbon — h- 
BN or SiC — Si02 nanoribbon — SiC one (Fig. la 
depicts a substrate with h-BN). The last two meth- 
ods produce a graphene sheet containing a nanoribbon 
without hydrogenation (as the nonhydrogenated one in 
Fig. lb) or a graphene strip without adsorbed CrOs 
molecules, respectively. Furthermore, the bandgap can 
be varied by using partially hydrogenated graphene 
(where some carbon atoms are not bonded to hydrogen 
atoms). Combinations of these methods can also be 
employed. Extensive experimental studies of graphene 
on substrates made of various materials, including rare- 
earth metals, have been reported recently [19-21]. It 
may be possible to open a bandgap via adsorption of 
other molecules on graphene or by using other materi- 
als as substrates. The use of gapped graphene to create 
potential barriers opens up additional possibilities for 
bandgap engineering in carbon-based materials [22]. 

We assume that both hetero junctions in the system 
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Fig. 1. Two configurations of the system under study: 
(a) graphene sheet on a substrate consisting of a Si02 
nanoribbon of width d inserted between h-BN 
nanoribbons; (b) a graphene sheet on a Si02 sub- 
strate containing a nonhydrogenated nanoribbon of 
width d, where open and closed circles are hydrogen 
atoms bonded to carbon atoms in different sublattices 
on opposite sides of the sheet, respectively. 



Fig. 2. First Brillouin zone of graphene, with linear 
energy spectrum at the corners (Dirac points). The 
reciprocal lattice vectors bi = (27r/3a, 2n j \JZa) 
and b2 = (27r/3a, —27r/V3a), where a = 1.42 A is 
the lattice spacing, combined with the dashed lines 
equivalent ly represent the Brillouin zone as a rhombus. 



combining a nanoribbon with gapped graphene sheets 
are type I junctions (e.g., see [23] for classification of 
junctions); i.e., the Dirac points in gapless graphene fall 
within the bandgaps of the adjacent gapped graphene 
sheets. This prevents spontaneous electron-hole pair 
creation, which would otherwise shunt graphene-based 
nanoelectronic devices such as FET. 

We believe that planar heterostructures made of 
gapless and gapped graphene are as prospective build- 
ing blocks in future carbon-based nanoelectronics. The 
use of only gapless graphene reduces the diverse op- 
portunities offered by bandgap engineering in gapped 
graphene. 

The paper is organized as follows. The preliminary 
remarks in Section 2 recall some aspects of the the- 
ory of quasi-relativistic fermions in graphene. Section 
3 introduces a generalized Dirac model for graphene. 
Section 4 describes the size quantization of originally 
massless carrier states in a planar graphene quantum 
well. Interface states are examined in Section 5. An ef- 
fective bandgap opening in the size-quantization spec- 
trum of the heterostructure in question makes it pos- 
sible to generate excitons by optical pumping. The ex- 
citon spectrum is calculated in Section 6. In Section 7, 
the effect of an applied uniform electric field on an ex- 
citon is analyzed. Experimental manifestations of the 
effects examined in this paper are discussed in Section 
8. The Conclusions section summarizes and discusses 
the main results obtained in this paper (Section 9). 

2. PRELIMINARY REMARKS 

In our previous study [24], we analyzed a planar 
heterostructure combining two narrow-gap semiconduc- 
tors with a graphene ribbon. We considered states with 



a definite parity A, which is an eigenvalue of the parity 
operator [25] 

P = *74A n - (1) 
Here, 274 is the inversion operator and 

A n = e-^ Sn = -zEn 

is the operator of rotation by tt about an n axis per- 
pendicular to the graphene plane. In standard repre- 
sentation, 

where a denotes Pauli matrices, and 
74-/3=( / _° / ) 

where / is the 2x2 unit matrix. It is clear that operator 
(1) is analogous to the parity operator i^n in quantum 
electrodynamics, where n = 711, 

7=(-°<xo), 75 = 71727374 = ^(?o). 

The eigenfunctions of this operator describe electron 
polarization states [26]. 

Charge carrier states in graphene can be described 
in terms of helicity define d as the eigenvalue of the 
operator h = cr • p/(2|p|). The projection of pseu- 
dospin on the direction of quasimomentum p indicates 
the valley in the Brillouin zone where electrons or holes 
belong (K or K' point in Fig. 2). Positive helicity cor- 
responds to electrons and holes with wavevectors near 
the K and K' points, respectively; negative helicity 
corresponds to electrons and holes with wavevectors 
near the K' and K points, respectively [27]. 

Massless states with opposite helicities are decou- 
pled [28]. In addition, charge carriers have chiral sym- 
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metry (helicity is conserved), and parity can be de- 
fined for both massless and massive carriers 1 . In other 
words, a higher symmetry of massless charge carriers 
implies the existence of an additional quantum num- 
ber: helicity. Whereas parity distinguishes only be- 
tween the valleys where carrier states belong (A = 1 
and -1 for states close to the K and K' points, re- 
spectively), helicity differs between a particle (electron) 
and an antiparticle (hole). However, chiral symmetry 
is broken for massive charge carriers (i.e., helicity is not 
a good quantum number any longer). Carrier states in 
a planar heterostructure combining gapless and gapped 
graphene should be characterized by parity. 

Recall that the Dirac equation describing massless 
carriers in graphene in terms of 4x4 matrices is de- 
rived by assuming that they are spinless and have two 
valley degrees of freedom [10]. When analysis is re- 
stricted to charge carriers in one valley, the Dirac equa- 
tion can be reduced to a 2x2 matrix representation 
by Weyl's equation for a massless fermion analogous 
to neutrino in two Euclidean dimensions. The car- 
rier energy spectrum with a pseudospin splitting in a 
planar heterostructure combining gapless and gapped 
graphene cannot be correctly analyzed in the 2x2 rep- 
resentation. For similar reasons, the representation of 
the Dirac algebra in terms of 2x2 matrices is not suf- 
ficient for describing the chiral symmetry breaking in 
quantum electrodynamics in two Euclidean dimensions 
[31]. 

Using the two-dimensional 4x4 Dirac equation to 
describe charge carriers in a graphene-based nanostruc- 
ture, we can study pseudospin effects following an ap- 
proach to narrow-gap semiconductor heterostructures 
based on the Dirac model [32]. This makes methods 
developed for solving problems in the spintronics of 
narrow-gap semiconductor heterostructures applicable 
to graphene-based ones [33-40]. 

As applied to a planar graphene-based heterostruc- 
ture, these methods can be used to manipulate valley 
occupation in graphene via control of pseudospin. A 
heterostructure of this kind can be used in valleytron- 
ics [41-47]. As a manifestation of pseudospin effect on 
charge carriers, we examine the pseudospin splitting of 
size-quantization spectrum (see Section 4 below). 

3. MODEL 

To describe size quantization in graphene-based het- 
erostructures, an equation containing a mass term should 
be written for the envelope wavefunction. A bandgap 
opening in the energy spectrum of graphene results 
from the lack of symmetry between the two triangular 
sublattices of its hexagonal lattice. The corresponding 
tight-binding Hamiltonian taking into account nearest- 
neighbor hopping has the form [48] 

h = -t Yl K( B + d ;)MB) + £ (BK(B + d;)] 



/A 



Fig. 3. Part of hexagonal lattice, with highlighted 
vectors df from a B sublattice atom to the three 
nearest-neighbor A sublattice atoms. 



B,cr 



di) - (B)MB)] , (2) 



where t « 2.8 eV is the nearest-neighbor hopping en- 
ergy; the sum runs over the position vectors B of all B 
sublattice atoms; the vectors (i = 1, 2, 3) pointing 
from a B sublattice atom to the three nearest-neighbor 
A sublattice atoms are expressed in terms of the lattice 
spacing a as 



d, = 



V3 



-a, a 

2 ' 2 



d 2 = 



d 3 = (-a, 0) 

(see Fig. 3); a =t, i is the (pseudo)spin index; a a (aj.) 
and b a (frj.) are the annihilation (creation) operators 
of A and B sublattice electrons, respectively; and the 
parameter A quantifies the on-site energy difference 
between the two sublattices (setting A = restores 
the symmetry between sublattices so that graphene be- 
comes gapless, whereas nonzero A equals the half-gap 
width in gapped graphene as shown below). 

Performing a Fourier transform, we change to the 
momentum representation 



a a (A) 



d 2 k 



a a (k)e 



>k-A 



mb) = / (£^ (k)eik ' B 



1 Massless states can be characterized by two quantum num- 
bers: helicity and sign of energy or helicity and eigenvalue of the 
operator 275 [29]. Parity is analogous to the eigenvalue of 275. 



where Qb means integration over the first Brillouin 
zone. 
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Hamiltonian (2) is rewritten as 
d 2 k 



(4(k) 6t(k)) 





(3) 



Conduction and valence band extrema lie at the cor- 
ners of the Brillouin zone. We use Hamiltonian (3) 
expanded around the K point with quasimomentum 

( 1 1 = (la' 3v5 a ) or aroun( ^ the point with c\2 = 

( 2tt 2tt V 

V 3a ' 3 V3a J ' 



H 



d 2 fc 

(27T) 



$t(^i(k). 



2 cr 



where integration is performed over small neighbor- 
hoods of the K and if' points. Near the corners, the 
Hamiltonian reduces to 



vpcr 



k- 





Aa z 



v F cr 




k 



■ Aa z 



(4) 



where = |at is the carrier Fermi velocity, a = 
(&x, &y) an d cr* = (cj Xl —(Jy) are Pauli matrices in the 
sublattice space, and ^(k) is the bispinor defined as 



*a(k) 



$S 1) (k) N 
$S 2) (k) y 



in terms of 



$(l,2) (k) = 



exp(— a z 



( «cr(qi,2 

V^(qi,2 



We consider the heterostructure combining gapped gra- 
phene regions 1 {x < 0) and 3 (x > d) with a graphene 
nanoribbon 2 (0 < x < d). The 2 axis is perpendicular 
to the graphene plane, the x axis is perpendicular to 
the hetero junction interfaces, and the y axis is parallel 
to the interfaces (see Fig. 1). 

Using Hamiltonian (4) , we write an equation for the 
envelope wavefunction in a planar heterostructure: 



[VFJ (t <8> CFxPx + T z <g) (JyPy) 

+ T <g> CT^Aj + T <g> CT (Vj - £)] 2/) = °- 



(5) 



Here, Aj = E g j/2 (j = 1,2,3) denotes half- width of 
bandgap (Ai ^ and A 3 ^ in regions 1 and 3, 
whereas A2 = region 2); the respective work func- 
tions V\ and Vs of regions 1 and 3 depend on the 
mid-gap energies relative to the Dirac points for the 
corresponding materials (we set V2 = to be specific, 
see Fig. 4); the 2x2 unit matrix cr acts in the sub- 
lattice space; the 2x2 unit matrix To and the matrix 
r z defined similar to the Pauli matrix <j z act in the 
valley space; ® is the Kronecker product symbol; and 
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J dx 



and p y 



ator components {h = 1) 



- are the momentum oper- 
Assuming that the carrier 
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Fig. 4. The quantum well under analysis. 

Fermi velocities may differ between the three regions, 
we denote those for gapped regions 1 and 3 by vpi and 
vfs and use vf2 ~ 10 8 cm/s for gapless graphene. 
Charge carriers move freely along the y axis: 

V(x,y) = ^(x)e ikyV . 

The wavefunction ty(x) is a bispinor: 



tf(a;) = 



where the spinors ^k(x) and ^k'(x) represent charge 
carriers in the K and K' valleys, respectively: 

V ' \$KB\x)j Y K ' \i) K 'B\x) 

In the present context, the parity operator is expressed 
as follows: 

P = T z <S)(T0. (6) 

Equation (5) is solved here in the parity basis. The 
eigenfunctions *ff\(x) of parity operator (6) are defined 
as follows: 

Py x (x) = \V x (x), 

^ +hK (x)^ 




Rewriting Eq. (5) as the 2x2 matrix equations 



*-l(2c) 



(7) 



-iw Fi CT x — + v Fj k y a v + XAj(j z + Vj J ^\k{x) ^ 
= E x iI>\k(x), 
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-IVpjC 



' dx 



V Fj ky(J y 



\AjCr z 



E\^\k> 0). 



(9) 

We see that setting Aj = and Vj = brings us back 
to the spinor wavefunctions describing chiral states near 
the K or K' point, where the operator h can be defined. 
However, chiral symmetry is broken when A / (see 
previous section) . Defining parity A as the eigenvalue of 
operator (6), we find that it indicates the valley where 
charge carriers belong: by virtue of (7), A = +1 for 
the states near the K point described by Eq. (8) and 
A = -1 for the states near the K' point described by 
Eq. (9). 

In both gapped and gapless graphene, the valleys 
transform into each other under time reversal. This is 
indicated by the opposite signs of the terms propor- 
tional to k y in Eqs. (8) and (9), since k y — » — k y under 
time reversal. It can be shown directly by using the 
time reversal operator T in explicit form that A — >• — A 
under T. Indeed, if k y is parallel to the line K — M — K' 
(see Fig. 2) and its origin is set at M, then K — >• K' 
and K' — >• K under T. 

Equations (8) and (9) are equivalently rewritten as 
the 2x2 matrix equation 



-iV F j(J X — + \V F jky(Jy 

: E\ip\ (x) . 



Aj(T z 



(10) 



Hereinafter, valley indices K and K' are omitted as 
unnecessary since A specifies the valley where charge 
carriers belong. 

We now discuss the boundary conditions at the in- 
terfaces between different graphene materials. At the 
outset, we note that they are easier to formulate than 
those used at the graphene-free-space interface in mod- 
els of edge states [6, 49]. To derive boundary conditions 
in the present model, we must find a relation between 
and ipx(-l) as I — >• in the neighbor-hood of x 
= 0, where I goes down to an atomic scale (condition 
at x = d is derived similarly). Multiplying Eq. (10) by 
ifj^(x) on the left, we integrate it over [— Z, I]. Since a 
is small, we neglect all terms except those containing 
a derivative with respect to x to obtain 2 

vl^H^M^H) = V>i +)t (04 + M +) (0, 

where ip^ ^ and are defined on the left- and right- 
hand sides of the boundary (at x < and x > 0, re- 
spectively). Representing these functions as 



exp (^iif^^j , 
we rewrite the equality above as 



(-) 



^- )f H) 



(+) 



^i +)t (0 



2 From the given equality, we have the continuity of the cur- 
rent component normal to the interface in the heterostructure 
plane. It is necessary condition. 



To formulate the boundary condition in final form, we 
assume that the difference between phases of i/)^ and 
near the interface is a multiple of 2tt: 

(/? (+) = (/? (-) + 27m, n e Z. 

As I goes to zero, we obtain the following wavefunct ion- 
matching condition [34, 35] 



where minus and plus signs refer to the materials on the 
left- and right-hand sides of the boundary, respectively. 

The solution to Eq. (10) is expressed as follows. 

1. At x < 0, 



'iPx(x) = C ( I e 



(ii) 



where 



.VFijh - Xky) 



VFiki = y Af - (E x - Vi) 2 + *' 2 
2. At < x < d, 

/ * \ 



V^^y. 



where 



1 j VFl 
VF2 



1 + 2 



c 



Xky 

~k^ 



(12) 



v F1 (k 1 - Xky) E x 



v F2 k 2 {Ex-V 1 + A 1 ) 



Q2 



VF2(k 2 + iXky) 
Ex ' 



E x = ±v F2 Jk 2 2 +kl 



with plus and minus corresponding to electrons and 
holes, respectively. 
3. At x > d, 



(13) 



where 

c = 



VF1 
VF3 



cos(k 2 d) 



Xky v F i(ki - Xk y )E x 
k 2 vf^Ex-Vx^t A x ) 



sm(k 2 d) 



VFzjkz + Xky) 
l E x -V 3 + A 3 ' 



v FS k 3 = ^A§ - (E X - V 3 ) 2 + v 2 F2 k 2 . 

The constant C is found by using the normalization 
condition for wavefunctions (11) — (13), 



oo 

J *{(x)* x (x)dx=l. 
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The carrier energy spectrum is determined by the 
dispersion relation 

tan(fc 2 d) = v F2 k 2 f(Xk y ;k 1: k 3: E x ) : (14) 

where 

f(Xk y ;k u k 3 ,E x ) = [v F1 (k! - Xk y )(E x -V 3 + A 3 ) 
+v F3 (k 3 + Xk y )(E x - Vl + Ai)] 
x [£ A (£ A - y x + Ai)(S A - F 3 + A 3 ) 
- v F2 v F3 \k y (k 3 + Xk y )(E x - Vi + Ai) 
+ v F iv F2 k y (ki - Xk y )(E x - V 3 + AA 3 ) 
-v F iv F3 (ki - Xk y )(k 3 + A/^)^]" 1 

is a function of &2 as well. Equation (14) must be 
solved for &2> and then the energy E\ is found. 

In the case of an asymmetric quantum well, the de- 
pendence of (14) on A gives rise to pseudospin splitting 
as the extrema of the dispersion curves shift away from 
Brillouin-zone corners. The dispersion relation predicts 
that E x (k y ) 7^ E_ x (k y ), and an energy splitting ap- 
pears near the conduction-band bottom at k y = fc* e : 

5E e s = lEt^J - E e +1 (k; e )\. 

A similar energy splitting appears near the valence- 
band top at k y = k* e : 

SE^ = \E^(k; h )-E^(k; h )\ 

Thus, a graphene nanoribbon becomes an indirect band- 
gap semiconductor analogous to silicon and germanium, 
where an electron-hole plasma can exist [50]. In the 
case of a symmetric quantum well (Ai = A3, V\ — V3, 
v F \ — v F3 ) band structure is invariant under parity 
and there is no pseudospin splitting [24]. 

4. SIZE QUANTIZATION 

Solving Eq. (14), we determine the size-quantized 
energies 

E Xb T (ky) = ±V F2 ^kl H (Xky)+k 2 y , 

where b T = 1,2,... labels electron (-) and hole (+) 
branches, respectively. The size-quantized energy spec- 
tra for symmetric and asymmetric quantum wells are 
shown schematically in Fig. 5. 

We now determine the carrier effective masses aris- 
ing because of size quantization in the graphene nanorib- 
bon in a planar heterostructure. Note that the effective 
masses are invariant under parity regardless of pseu- 
dospin splitting. Hereinafter, we omit indices re- 
stricting ourselves to a particular branch of the electron 
spectrum and a particular branch of the hole spectrum. 

We write the dispersion law for electrons near an 
extremum at Afc* as 

y c 

^^E^^(k y -xk; e )\ 

1 A/ K 20e + K ye 
771 = 

e v F2 1 + k' 2 \ e + k 20e k f 2 0e ' 




Fig. 5. Energy spectra: (a) symmetric quantum well 
(no pseudospin splitting), with matching branches for 
A = +1 and A = -1 (E x (k y ) = E_ x (k y ) = E x (-k y )); 
(b) asymmetric quantum well, with pseudospin split- 
ting manifested by the "spread-out" in quasimomen- 
tum between the extrema at k* e for electrons, shown 
for b- = 1, and at k* h for holes, shown for b + = 1 

(E X (ky) ± E- X (ky)). 

where the respective values k 2 o e , k 2 Q e , k 2 f Q e of k 2e (k y ) 
and its first and second derivatives at k y = A/c* e are 

independent of A; Eq = v F2 ^k 2 0e + k* 2 is the energy 
at the extremum. 

Analogous expressions are obtained for hole ener- 
gies: 

E h x~E% + ^(k y -\k; h ) 2 , 

l k 2 +k *2 ( 16 ) 
1 y K 20h + K ye 

where the respective values k20h, k' 20h , k'^ of k<ih(ky) 
and its first and second derivatives at k y = — Xk* h ; 

Ek = -v F2 ^/k* 0h + kf h . 

To estimate characteristic values, we consider the 
planar heterostructure combining a gapless nanorib- 
bon with gapped graphene sheets with Ai = 0.75 eV, 
v F i = l.lv F2 , A3 = 1 eV, and v F3 = 1.2v F2 . The 
nanoribbon width is d — 2.46 nm (ten hexagonal cells). 
Since the unknown values of V± and V 3 can be found by 
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(a) 













(b) 






m e /m o 




Fig. 6. Electron and hole effective masses in the Fig. 7. Electron and hole extremum energies and in 
graphene nanoribbon (in units of free-electron mass the size-quantization spectra as functions of V\ for V3 
mo) as functions of V\ for V3 = (a) and as functions = (a) and as functions of V3 for V\ = (b). The 



of V 3 for V ± = (b). 



effective bandgap E e J * = Eq 
insignificantly. 



\ E o\ 



629 meV varies 



comparing our results with experimental data, we seek 
the dependence of energy spectrum parameters on V\ 
and V 3 . Note that |Vi| < Ai and \V 3 \ < A 3 to ensure 
that the heterostructure is type I. 

Figures 6-9 show the results of numerical calcula- 
tions of electron and hole effective masses in the graphe- 
ne nanoribbon, extremum energies, k* e and k* h val- 
ues, and pseudospin splitting SEp h plotted versus work 
function for one of the gapped graphene sheets given 
that the work function for the other is zero. 

It is clear from Fig. 9 that the pseudospin split- 
ting energy may amount to approximately 10 meV. To 
obtain a larger pseudospin splitting, the quantum well 
must be more asymmetric. Both V± and V3 can be var- 
ied by shifting the valley energies in gapped graphene 
under applied stress, with potential barriers playing 
the role of bandgaps in the gapped graphene sheets. 
An analogous effect is achieved by applying an electric 
field on the order of 10 6 V/cm perpendicular to the 
interfaces in the graphene plane [40]. 

As expected, the energy spectrum is symmetric un- 
der the change E — > -E when V\ and V3 = 0; i.e., the 
electron and hole spectra have equal effective masses, 
extremum energies, extremum positions, and pseudo- 
spin splitting energies. The electron and hole effec- 
tive masses in graphene are smaller than those in the 
gapped graphene sheets adjoining the gapless graphene 
nanoribbon (m* = Ai/v^ « O.llrao and = A 3 /vp 3 
« 0.15ra ). 



5. INTERFACE STATES 

The existence of surface states was predicted by 
Tamm in [52]. Extensive studies have been conducted 
of Tamm states in various systems including semicon- 
ductor superlattices [53, 54]. We consider interface 
states of a new type that arise in a narrow quasimo- 
mentum interval from the crossing of dispersion curves 
and are analogous to those in narrow-gap semiconduc- 
tor heterostructures [55, 56]. In the planar graphene- 
based heterostructure examined here, these states are 
localized near the hetero junction interfaces between 
the nanoribbon and the gapped graphene sheets. Inter- 
face states can exist not only in quantum wells but also 
in quantum barriers [10]. Note that interface states 
arise as well from the crossing of dispersion curves in a 
single hetero j unction between different graphene ma- 
terials [57]. 

The wave function describing an interface electronic 
state is expressed as follows. 
1. At x < 0, 

where 

vfiKi = yjAl-iEx-Vtf + vU*. 
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^—■i= +1 

iX P — i 


A = 





100 

V , meV 



Fig. 8. Extremum points of size-quantization bran- 
ches for electrons (fc* = fc* e ) and holes (fc* = k* h ) as 
functions Vi for Vs = (a) and as functions of V3 
for Vi = (b). Inserts show the relative positions of 
dispersion curves for Vi = Vs = (a) and at fc* e = 
z,. 

for simplicity. 



yh (b); if and K' points are set at the same position 



Fig. 9. Pseudospin splitting in electron and hole spec- 
tra, 8 El and 5E^ , as functions V\ for V3 = (a) and 
as functions of Vs for Vi = (b). Vanishing SEp h 
corresponds to vanishing k* eh in Fig. 8, as shown in 
inserts to (a). Inserts to (b) show positions of disper- 
sion curves when /c* e and k* h coincide. 



2. At < x < d, 



C 



where 



K± 



1 I vfi 

2VVF2 



1 ± 



Xko, 



± 



v F i(ki - ^k y )E x 



K2 V F2 K2(E X -V 1 +A 1 ) 



.VF2^2 + A&J ^ .V F 2^2 - Xk y ) 

Q2 = I ^ ~ > ^2 = 



E\ 



E\ = ±v F2 \/kl - k%, 



(17) 



with plus and minus corresponding to electrons and 
holes, respectively. 
3. At x > d 



where 

c = 



VF1 
VF3 



Gh(K 2 d) + ' ^ 



^ F i(^i - \k y )E x 

K2 ' V F 2K2(E X ~ Vi + Ai) 



sh(n2d) 



Q3 



E\ — Vq 



^F3&3 



(£7a - ^) 2 



^F3 



The expression for energy in (17) implies that an 
interface state exists only if 3 

\K 2 \ < \ky\. 

We obtain the dispersion relation 

tanh(« 2 d) = v F 2K2f(^k y ] «i, « 3 , E\). (18) 

which is similar to (14) up to the substitutions k\ — >• 
^1, k.2 — >• ^2, and /C3 — » ^3. 

We calculated numerically the energies of interface 
states as functions of quasimomentum component k y 
for heterostructures with Ai = 0.75 eV, A3 = 1 eV, Vs 
= 0, vfi = v F s = 0.7^F2, and Vi = and 100 meV. The 
results shown in Figs. 10 and 11 demonstrate that only 
electron interface states with A = +1 and hole interface 
states with A = -1 occur. 

When Vi = 0, the allowed quasimomenta for hole 
interface states (A = -1) are similar to those for elec- 
tron states, but the hole and electron energies have 
opposite signs; i.e., the spectrum is symmetric under 
the change E\ — » —E\. When V\ = 100 meV, the sym- 
metry is broken and hole interface states exist only at 
negative quasimomenta. 



3 The zero mode corresponding to |«2| = \k y \ (with E\ = 0) 
is irrelevant here because i/j\(x) = 0. 
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dk 



-2.3 -2.2 




2.9 -2.6 -2.7 -2.6 




2.55 2.6 



Fig. 10. Energy of electron interface states with A = 
+1 in the heterostructure with Ai = 0.75 eV, A3 = 
1 eV, vpi = VF3 = 0.7vp2, and V\ = V3 = (electron 
interface states with A = -1 do not exist). Allowed 
quasimomenta lie in two intervals: (a) negative k y and 
(b) positive k y . 

6. EXCITON IN A PLANAR GRAPHENE 
QUANTUM WELL 

In gapless graphene, the carrier effective mass is 
zero and excitons do not exist. The existence of ex- 
citons in gapless graphene would lead to excitonic in- 
stability and excitonic insulator transition to a gapped 
state [58, 59]. 

The energy gap arising in a graphene nanoribbon as 
described in Section 4 makes it possible to generate ex- 
citons by optical excitation or electron-hole injection. 
Quantum-well excitons strongly affect optical proper- 
ties of the system considered here. 

Excitons in similar quasi-one-dimensional carbon- 
based systems (semiconducting single- and multi- walled 
nanotubes) have been studied theoretically in [60] . The 
exciton spectrum is calculated here for a planar graphene 
quantum well by using the model applied to quantum 
wires in [61]. This model yields simple analytical ex- 
pressions for exciton binding energy. 

Since formulas (15) and (16) are obtained in the 
nonrelativistic limit, the two-particle exciton wave func- 
tion depending on the electron and hole coordinates y_ 
and y+ in a sufficiently narrow nanoribbon must obey 
the ID Schrodinger equation with Coulomb potential: 



> 

p 




Fig. 11. Same as in Fig. 10, but with V\ = 100 
meV; (c) energy of hole states with A = +1 for allowed 
quasimomenta in a single interval of negative k y . 



where E' — E — E e /f and e 2 = e 2 /n eff . The effective 
dielectric constant of graphene, n e ff = (e + e')/2, may 
vary widely with the dielectric constants e and e' of 
the media in contact with graphene, such as free- space 
permittivity and substrate dielectric constant [62, 63]. 

The electron-hole Coulomb interaction in a ID gra- 
phene nanoribbon is three-dimensional, but the prob- 
lem can be reduced to one dimension (electron and hole 
y positions) for sufficiently narrow nanoribbons. 

Rewriting Eq. (19) in terms of electron-hole sepa- 
ration y = y_ — y+ and center-of-mass coordinate 



Y 



KV- + ™* h y+ 
mi + ml 



1 d 2 



1 d 2 



Z 2 



2m* dy 2 _ 2m^ dy\ 
E'(/>(y-,y+), 



\y- -y+\ 



<f>(y-:y+) 



and introducing the function 

<j>(y-,y + )=Myy KY , 
where K is the total exciton momentum, we obtain 



(19) 



2/x* dy 2 \y 



ip n (y) = E n ip n (y), 



(20) 
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where /i* = m*m^/ (m* H-raj£) is the reduced mass and 
E n is the energy of the nth exciton level (n — 0, 1, 2, ... 
is the principal quantum number). The total exciton 
energy E' is obtained by adding the total kinetic energy 
of the electron-hole pair to E n : 



E' — E n 



K 2 



2(m*+m*) 



To find the solution at y > 0, we substitute ^ n {y) 
represented as 



i>n(y) = Bn exp (-y/a n ) F n 



2y 



into Eq. (20) and obtain the confluent hypergeometric 
differential equation 

^-^ + V F n = 0, (21) 

where £ = |^ and r\ = /i*e 2 a n . We also have 

1 



EL 



(22) 



Equation (21) with r] = n is solved by the associated 
Laguerre polynomial 



1 jn 

n! a£ n v 



»-i e -q 



and the wavefunction is expressed as 

^n(2/) = 5 n exp(-?//a n )L^ 1 ( — ) . 

Analogously, we find the solution to Eq. (20) at 
y < 0: 



^n(y) = ±B n exp(y/a n )L n 1 



2y 



where and "-" are taken for n = and n / 0, 
respectively, and the continuity of ip n (y) and its first 
derivative ijj' n (y) are used as boundary conditions. Since 
VVi(0) = and ^'(0) 7^ for n ^ 0, the excited-state 
wavefunction ip n (y) is odd (otherwise, it would be dis- 
continuous at the origin), whereas the ground-state 
wavefunction is even. 

The normalization condition 



00 

/ 



\i> n (y)\ 2 dy = l 



is used to determine the coefficient B n in the expression 
for ij) n (y)\ 



00 



-1/2 



and B n = \/y/2a n for n = 1, 2, . . . and Bq = 1/y/ao 
for n = 0. 



Here, 



(23) 



(n = 1, 2 . . .) is the Bohr radius of an exciton in the 
nth excited state. Combining (22) with (23), we find 
the exciton energy spectrum: 



/re* 



(24) 



The ID ground-state (n = 0) Coulomb energy ex- 
hibits a logarithmic divergence at short distances [64]. 
Therefore, the lateral spread of the exciton wave func- 
tion (along the x axis) due to the three-dimensional 
nature of Coulomb interaction should be taken into 
account by introducing a cutoff parameter do < d. Av- 
eraging the kinetic energy operator 



2/i* dy 2 



and the potential 



V(y) 



v 



'(M-do) 



over ground-state wave functions 

My) = -Le-w/ ao , 



(25) 



where the ground-state Bohr radius ao plays the role 
of a variational parameter, we express the ground-state 
exciton energy as 



1 2e 2 a 

Eo = - 2 In — 

2/i*aQ ao a 



(26) 



Minimizing (26) with respect to ao, we obtain an equa- 
tion for ao: 

(27) 



a 



ai 



2(lnf -I)' 
To logarithmic accuracy, when 



we find the relations 



E = 4£i In 



2 ai 
d' 



a 



ai 



2 In 



(28) 

(29) 
(30) 



Using (27), we easily obtain the next-order correction 
to E : 

^ = -8^ In I In (2 In |). 

We now examine the applicability of the formulas de- 
rived here. The semiconducting state induced in a 
graphene nanoribbon is stable with respect to spon- 
taneous electron-hole pair creation (excitonic insula- 
tor transition) only if the exciton binding energy \Eq\ 
is smaller than the effective bandgap in the graphene 
nanoribbon, 

\E \ < Efl. 
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Furthermore, the quantum well width d must be much 
smaller than the exciton Bohr radius ai, 

d <C a\. 

Logarithmically accurate formula (29) is correct only if 
condition (28) holds. However, the asymmetric quan- 
tum well analyzed here to examine pseudospin effects 
may not admit even a single size-quantization level if 
the graphene nanoribbon width d is too narrow. As 
d decreases, the effective bandgap increases, approach- 
ing A + + A_, where A± = min{Ai ± V u A 3 ± V 3 } 
(with plus and minus corresponding to electrons and 
holes, respectively). When a certain d c is reached, the 
size-quantization levels are pushed into the continuum. 
This imposes a lower limit on d: 

d > d c , 

where d c can be estimated as 

. 7TV F 2 



c A + + A_ 

As d increases, condition (28) is violated. In this case, 
a more accurate variational calculation should be per- 
formed using the modified three-dimensional Coulomb 
potential 



V(y) 



g 2 



Vv~ 2 



di 



where do is a cutoff parameter. We average the Hamil- 
tonian with potential V(y) over trial functions (25) to 
obtain 

1 2^ 
£0 = - — /(p), (31) 



2/i*aQ 



a 



where I{p) = | [Ho(p) — Yo (/?)], H v (p) is a Struve func- 
tion, Y u (p) is a Bessel function of the second kind, and 
p = 2do/ao (y is here and 1 below). 

Minimizing (31) with respect to ao, we obtain an 
equation for ao: 



2a ( , 4d ( , 



1, 



(32) 



where J{p) = 1 - § [H^p) - Y 1 (p)}. 

Figure 12 shows the numerical results obtained by 
using both methods to calculate Eo(d) for the het- 
erostructure described in Section 4, with do oc d ad- 
justed to match the curves at small d. Discrepancy at 
large d increases as ln(ai/d) approaches unity. 

7. ELECTRIC FIELD EFFECT 
ON EXCITON LEVELS 

Interaction between an exciton and an external elec- 
trostatic field £ is described by the operator 

Hi -dS = \e\(£ x x + £ y y) 

where x = \x- — x+\ and y = \y_ — y+\ are the elec- 
tron-hole relative position vector components and d 
is dipole moment. The electric field is supposed to be 
weak enough to ensure that the energy level shift is not 




Fig. 12. Exciton ground-state energy calculated by 
formula (29) (dashed line) and by formula (31) after 
Eq. (32) is solved numerically for ao (solid line); do = 
0.22d. 

only smaller than the spacing between size-quantization 
levels but also smaller than the spacing between exci- 
ton levels. These conditions can be written as 

d <C ai <C as, 

where ag = (/i*|e|£) -1 / 3 is the electric length. 

We consider two cases: (1) the electric field is ap- 
plied parallel to the x axis and perpendicular to the 
nanoribbon edges in the graphene plane; (2) the elec- 
tric field is applied along the y axis, parallel to the 
nanoribbon edges. 

In the former case, the energy shift varies linearly 
with the difference between the average x components 
of the electron and hole position vectors: 



j- w — 



|e|£«a:_>A -<*+>*)> 



(33) 



where average x components are calculated by using 
electron and hole single-particle wave functions, gen- 
erally depending on the electron and hole eigenvalues 
A and A' of the operator P, respectively. Exciton en- 
ergy shift (33) is independent of the principal quan- 
tum number n. It may vary with A and A 7 , resulting in 
different exciton binding energies (more precisely, the 
binding energy of an electron-hole pair with A = ±1 
and A' = ±1 may have four different values). 

In the latter case, the first-order electric field-induced 
correction is zero 4 , 



i\£(y)n = 0, 



(34) 



because the integral of y\^ n (y)\ 2 with respect to y 
vanishes. To evaluate the second-order electric field- 
induced correction, we make use of the Dalgar no-Lewis 
perturbation theory [65]. Defining a Hermitian opera- 
tor such that 



[F, H ]\n)=Hi\n) 



(35) 



where 



2/i* dy 2 



4 Note that E 



(i) 



in the former case if the electron and 



hole spectra transform into each other under field inversion. 
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is the zeroth-order Hamiltonian, \n) = ipn(y) is the 
zeroth-order wave function of the nth exciton level, and 
Hi = \e\£y, we obtain 



r(2) _ 



(n|^F|n)-(n|^|n)(n|F|n). (36) 



In the case in question, the second term in this formula 
vanishes by virtue of (34). 
Rewriting Eq. (35) as 



d 2 F d^ n dF 



dy 2 



dy dy 



(37) 



we find 



F(y) = 2 M * 



dy' 



\Mv')\ 2 

— oo 

y' 

x J dy"r n (y")HnPn(y")- 



(38) 



Combining (36) with (38), we have the exciton gro- 
und-state energy shift 



(2) 



128 In 4 a 



(39) 



which is very small compared to Eo given by (29) be- 
cause of the fourth power of a logarithm in the denom- 
inator and a small numerical factor. 

For comparison, we write out the energy correction 
to the first excited exciton state: 



(2) 



;(31- 6 7 )a?£ 2 , 



where 7 = 0.577. . . is Euler's constant. 

By analogy with layered heterostructures [66], the 
ionizing (exciton-breaking) field strength 6 C is estimated 
as 

\E,\ 



£c 



(40) 



m(\y\w 

where (\y\)o = ao/2 is the average electron-hole sep- 
aration for the ground-state exciton. To logarithmic 
accuracy, it follows that 



Sc 



*2 i~ 

fi \e 



5 ln 3 %. 
a 



(41) 



To get the order of magnitude of £ c , consider the 
quantum well discussed in Section 4. Setting m* = 
mj* ~ 0.0056mo, the Si02 substrate dielectric constant 
n e ff « 5, d = 2.46 nm, and a\ ~ 81 nm, we use for- 
mula (41) to obtain £ c = 9 kV/cm. 

8. POSSIBLE EXPERIMENTS 
ON THE HETEROSTRUCTURE 

Pseudospin splitting can be observed by means of 
Raman spectroscopy. The D' peak of interest for the 
present study (alternatively called 2D peak to empha- 
size that it is due to a two-phonon-assisted process) 
is located at 2700 cm 1 [67]. It arises from intervalley 




Fig. 13. Possible double resonant Raman processes 
involving electron scattering between valleys. To sim- 
plify presentation, analogous processes involving hole 
scattering between valleys are not shown. 



scattering involving phonons with wavenumbers q> K, 
where K = 4ir/3y/3a « 1.7 x 10 8 cm 1 is the spacing 
between adjacent K and K r points. One process of this 
kind is indicated a,sA^B^C^D^Am Fig. 13. 

Pseudospin splitting enables intervalley scattering 
involving phonons with q' w q =p Ak (with plus for 
electrons and minus for holes), where Ak = 2/c* e and 
Ak = 2k* h in electron and hole scattering, respectively. 
These processes contribute to a peak blueshifted from 



D' by Auj^ and a peak redshifted from D' by Auj^ r ^ , 
giving rise to a doublet structure of the D' peak. 

An estimate for Auor can be obtained by using op- 
tical phonon dispersion uj p h(q)- The Raman shift is 
twice the optical phonon frequency: 

6uj R (q) = 2u ph (q). 

The change in the Raman shift caused by pseudospin 
splitting is 



°R 



\5uj r (AKt^)-5uj r (AK)\. 



which amounts to Aujr w 24 cm -1 for characteristic 
values of the heterostructure parameters. This value 
essentially exceeds the Raman spectral resolution of 
1 cm 1 [68] and compares to the D' peak width for 
gapless graphene, To = 30 cm -1 [69, 70]. 

Note that a blue shift of the D' peak has also been 
observed in the Raman spectrum of epitaxial graphene 
on a SiC substrate [68] . This effect is attributed to the 
strain induced by the substrate in quasi-free graphene 
since the SiC lattice constant exceeds substantially that 
of graphene. 

Raman scattering contributions from gapped gra- 
phene sheets can be avoided either by using a laser 
beam whose width is smaller than that of the gapless 
graphene nanoribbon (d < 10 nm) or by pumping at 
a frequency lj such that the beam cannot be absorbed 
by gapped graphene materials, 



E eff 



2u p h < uo < rain{2Ai, 2A3} 



The positions of the luminescence lines correspond- 
ing to exciton levels can be determined from optical 
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experiments and compared to theoretical predictions. 
The splitting of exciton lines in an electric field is eval- 
uated by using formulas (33) and (39), respectively. 

Let us now discuss the Landau levels in a planar 
graphene heterostructure induced by a magnetic field 
perpendicular to the graphene plane. The magnetic 
field is supposed to be relatively weak, 



an = \/c/\e\H > d, 

so that size quantization plays a dominant role. 

As in analyses of layered narrow-gap heterostruc- 
tures [38, 71], the quasimomenta of free- moving carri- 
ers in size-quantized spectra (15) and (16) are replaced 
by the quantized magnetic momentum 



kHn 



2n 



a H 



where n = 0, 1, 2 . . . labels Landau levels. (While 
each size-quantization level splits into a respective set 
of Landau levels, the present analysis is restricted to 
the lowest size-quantization level.) 
Then, we have the carrier energies 



rph 
E n \ 



(Vn- 



where 



= ^-^(v^ + AvW) 2 , 

m* H c 1S the cyclotron frequency and 

n H e,h = a 2 H k;lJ2. 

The minimum energy of electrons (near K point, 
with A = +1) or holes (near K' point, with A = -1) 
may match a Landau level with n / 0. The zeroth 
Landau level is lowest only if riHe,h < 1/4. When 

N n -i < n He ,h < N n , 



where N n = [y/n + 1 + y/n) /4, the minimum energy 
(ground state) matches the nth Landau level. 

Such an unusual ground state can be observed only 
at temperatures sufficiently low that the excitation en- 
ergy from the lowest Landau level to the next one 
satisfies the condition AE e ^ h > T. It can easily be 



e,h 
H 



shown that AE\ 
ima at n Hmin = N n (n 

riHmax = (y/n + 1 + V 7 



H 

as a function of nn e ,h has min- 
0, 1, 2 . . .) and maxima at 
n - I) 2 / 4 (n = 1, 2 ...): 



AE e jj 



,h 



- 1 - ^n Hejh ) 2 - [y/n - y/n H e,hY 



N n -\ < TlHe,h < riHmaxi 
< TlHmax 



(y/n + l - ^/n H e,h) 2 ~ {Vn - ^n He ,h) 2 
< n He ,h < N n . 



A degenerate ground state corresponds to zero excita- 
tion energy. The excitation energy AE^ h is plotted as 
a function of riHe,h m Fig. 14. 

Finally, interface states can manifest themselves in 
the I-V curve of the planar heterostructure carrying 
a current parallel to the gapless graphene nanoribbon. 




Fig. 14. Excitation energy AE e ^ h vs. riHe,h- 
The insert shows positions of Landau levels with 

Ni < n He ,h < N 2 . 



An increase in applied electric field may cause charge 
carriers to "drop" into interface states (preferable ener- 
gy-wise), giving rise to a region of negative differential 
conductivity in the I-V curve. 

8. CONCLUSIONS 

We have analyzed the characteristics of planar gra- 
phene nanostructures (quantum wells). On the one 
hand, they retain the unique properties of infinite gra- 
phene sheets. On the other hand, bandgap opening 
makes them important building blocks in carbon-based 
nanoelectronics, which can be used to control electron 
motion. 

Parameters of graphene quantum wells can easily be 
manipulated by varying the gapless nanoribbon width 
or the potential barriers in the adjoining gapped gra- 
phene sheets. We predict pseudospin splitting to occur 
in asymmetric graphene quantum wells and interface 
states to arise from the crossing of dispersion curves of 
gapless and gapped graphene materials. 

We have performed calculations of optical proper- 
ties of planar graphene nanostructures and suggested 
possible experiments to study the effects in question. 

Analysis of pseudospin (valley) characteristics in 
the heterostructure is simplified by using an effective 
Hamiltonian having a pseudospin-split energy spectrum. 
Note that an analogous spectrum was discussed in [72, 
73]. Therefore, the effective Hamiltonian must contain 
a Bychkov-Rashba spin-orbit coupling Hr = ctR [crp] v 
responsible for pseudospin splitting, where olr is a Rash- 
ba parameter and v is the unit vector normal to the 
hetero junction interface in the graphene plane. In the 
coordinate system used in this paper, effective electron 
and hole Hamiltonians can be written as 



eff 



Py 
2m* 



Xa Re a z p y + A e 
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and 

p 2 

H eff = + Xa Rh°zPy + Aft, 

where = /c* e /m* and = ky h /m^ re the re- 
spective electron and hole Rashba parameters (which 
are different when the symmetry under field inversion 
is broken) and the parameters A e ^ are expressed in 
terms of extremum energies and Rashba parameters as 

A e = El + ^m* e a 2 Re , 
&h = -E% + ^m* h a 2 Rh . 

For characteristic values of the heterostructure param- 
eters, our estimates predict aR e ^ ~ 5 • 10 -8 eV-cm, 
which is higher by two orders of magnitude than for 
GaAs [72] (mainly because of small effective mass). 
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